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Abstract 



We discuss the Schrodinger functional in lattice QCD with staggered fermions 
including its order 0{a) boundary counterterms. We relate it, in the classical 
, continuum limit, to the Schrodinger functional as obtained in the same limit with 

Wilson fermions. We compute the strong coupling constant defined via the Schrodinger 
' functional with staggered fermions at one loop and show that it agrees with the 

■ continuum running coupling constant in the Schrodinger functional formalism. 

0^ 



1 Introduction 

From the pioneering work of the ALPHA collaboration is has become clear that 



^ ■ the Schrodinger functional in lattice QCD is a useful setup and tool for nonperturbative 
computations in QCD. Through the boundary conditions in the Schrodinger functional one 
can introduce, in a gauge invariant way, non-vanishing background fields and use those 
^ ! to probe the model. The response to a constant chromoelectric background field, for 
5^ \ example, allows the nonperturbative computation of a well-defined renormalized running 
strong coupling constant with the scale given by the spatial size of the system [jl]-^,^. 
The Schrodinger functional has also proved useful for the nonperturbative computation of 
the "clover" coefficient for 0{a) improved Wilson fermions and the gauge independent 
computation of current renormalization constants [0]. 

In this paper we shall discuss the formulation of the Schrodinger functional with 
staggered fermions. The basic steps have already been done by Miyazaki and Kikukawa 
We shall extend their work in several directions: discuss the form of 0{a) boundary 
counterterms, compute the one-loop fermionic contribution to the pure gauge 0{a) boundary 
counterterm tT{Fk4:FkA), and most importantly, discuss the relation between the lattice 
Schrodinger functional with staggered fermions and the one with Wilson fermions in 
the continuum limit. In particular we will compute the contribution from the staggered 
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fermions to the running coupling constant defined through the Schrodinger functional at 
one loop and show that it agrees with the computations using other regularization schemes. 
As a by-product of this computation we obtain, we believe for the first time, an explicit 
calculation, at one loop, of lattice artefacts from staggered fermions in a physical observable. 

In the next section we will briefly review the formulation of the Schrodinger functional 
with staggered fermions. We shall consider massless fermions, and discuss inclusion of a 
mass term later in section |[ In section |^ we will discuss the Schrodinger functional in 
the continuum limit in terms of four flavors of Dirac fermions. In section ^ we will compute 
the one-loop contribution from the staggered fermions to the Schrodinger functional running 
coupling for massless fermions. In section |^ we discuss inclusion of a mass term and compute 
the one-loop contribution to the coupling from massive fermions. Section ^ contains some 
concluding remarks. The symmetry properties of staggered fermions are briefly reviewed in 
appendix A, and a "two-time-slice" transfer matrix for staggered fermions is sketched in 
appendix B. 



2 The Schrodinger functional with staggered fermions 

The Schrodinger functional describes the evolution of a state at (Euclidean) time t = to 
another state at time t = T. Using the transfer matrix it can be written as a path integral 
with fixed boundary conditions at time t = and T. For fermions, described classically by 
a first order differential equation, actually only half the degrees of freedom can be specified 
at each boundary |^,10|. The transfer matrix for staggered fermions has been worked out 
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TT[| . It is non positive, but its square is positive. This leads to a doubling of degrees of 
freedom and it is better to think of this square as a transition amplitude from two time-slices 
to two time-slices (see appendix B for a "two-time-slice" transfer matrix that corresponds 
explicitly to the square of the transfer matrix derived in []TT|). As shown in this leads to 
the fact that for staggered fermions all degrees of freedom can be fixed at both boundaries. 
The Schrodinger functional can thus be represented as the path integral 

T— 1 

Z[WXX;W',C,C] = f [DU] I nil mx,x,)dxix,x,)]eM-SG - Ssf}. (2.1) 

X X4 = l 

Here W and W represent the boundary values of the gauge fields, Sq is the pure gauge 
action and [DU] the Haar measure over gauge fields, both appropriate for the Schrodinger 
functional formulation as in [ffl]. The fermionic part of the action is given by H] 



T-l 

^SF = Y.mi 2<^)x{x) [U,{x)x{x + /x) - Ul{x - ^i)x{x - /i)] + + 5^ (2.2) 



X X4 = l ^ 



with rj^{x) = (— 1)S''<M^"' the usual staggered phase factors. At the boundaries the fields 
take on their boundary values 

x(x,0) = C(f), x{x,0) = ax), x{x,T)=C{x), x{x,T)=C{x), (2.3) 
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with (, (, (' and (' independent complex Grassmann fields. The additional boundary terms 
of the action are 

3 ^ 

X k=l 

+ 5^^r/4(f,0)C(f)x(x,l). (2.4) 



2 

X 



and 



3 . 
X k=l 

- 5^^r74(x,T)C'(x)x(x,T-l). (2.5) 

X 

where Wk and Wj^ are the gauge fields at the boundaries. 

Note that for staggered fermions the total number of time-slices has to be even. 
Therefore, in the labeling adopted in this paper for the lattice sites, which agrees with 
that of the time extent T has to be odd. In the spatial direction we take the lattice to 
be of size L (even!) and impose the generalized periodic boundary conditions |^ 

Xix + Lk) = e^'^xix), Xix + Lk) = xix)e-''K (2.6) 

These boundary conditions are easily implemented by transforming to periodic fermion fields 
in a constant abelian background field = 1, Uk = q!'^^/^ . 



3 The Schrodinger functional in terms of four-component spinors 

We construct four- component spinors from the one-component Grassmann fields x ^"^^ X 
the standard way, following Kluberg-Stern et al. []T2[. We discuss here only free staggered 
fermions. Including the coupling to gauge fields is quite straightforward, though notationally 
somewhat cumbersome p!^JT3[] . Since we will be interested in the classical continuum limit, 
we shall display the lattice spacing explicitly (except in the labelling of the lattice sites). We 
divide the lattice into 2^ hypercubes in which the four-component spinors reside.Q Thus we 
set a; = 2y + ^ with ^^ = 0,1 and define 

x{2y + = x^y) , x(2y + = x^{y)- (s.i) 

We also introduce 

T^ = lhflM'- (3.2) 



^This division of the lattice into 2^ hypercubes requires that in each direction we have an even number 
of sites. With our labeUng convention in the time direction T is odd and < 2:4 < T. 
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We use hermitian Euclidean Dirac matrices. The matrices satisfy the orthogonahty and 
completeness relations 

tr (^r|rg,) = 4%,, J2 (r?)"" r/'' = as'^'S''^ (3.3) 

The four flavors of four-component Dirac spinors are then constructed as 

r^iy) = lJ2^rxdy)^ r"iy) = lJ2xdy){^y"- (3.4) 

where greek superscripts denote spin indices and roman superscripts flavor indices. These 
relations can be inverted: 

x^(y) = 25^(r|)'^"^"'^(y), xe(y) = 2 5^^'^"(2/)rr- (3.5) 

aa aa 

It is useful to introduce some more notation. Let A = (g) for some gamma matrix 
acting in Dirac space and some gamma matrix Tp acting in flavor space. Then 

A1.A2 = (r^.r^J ® (r^.r^J^ 
(A.r^)"'^ = (r^r^r^)"'' (3.6) 
(rJ.A)'^" = (ri^rjis)"". 

In the second and third equations above, just as in the definition of the ■?/'- fields, (|3.4| ), is 
a "mixed" matrix with the first index a spinor index, and the second one a flavor index. We 
also define 

i = l®l, f^ = 7^®l, = 75 ® (757m)^ (3.7) 

and the projectors 

1 r 1~ 



PI,^^ = - [1 ® 1 + (7^75) ® (757^f ] = -f ^ (f ^ + f 

Pi'^ = ^ [1 ® 1 - (7m75) ® (757,)T = , (f, - fj) . (3.8) 
These projectors are useful since 

^0^"^ • = 5c.,or5, rt . p(^) = 5^^,orJ, p^^^ ■r^ = 6^^,,r^, rl ■ pi^^ = 6^^,,rl. (3.9) 

Thus they project onto the one-component fields with = or 1. 

Inserting the relations (|3.5|) into the free staggered action, it can be written as 



Sp = (2a)^^^(y) \t^D^^{y)+tlaA^ij{y)] (3.10) 
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where 

^f^fiy) = -^[fiy + f^)-2fiy) + fiy-fi)]. (3.11) 
The expression for the action can be compactified further by defining as 

V^^Piy) = (1 ® l)D^^ljiy) + (7^75) ® i^.^^faA^^Piy). (3.12) 

To deal with the boundaries in the Schrodinger functional we note that the projectors 
Po,i = Pq'^i project onto the boundary fields 

PomO) = p{y), mO)Po = p{y), PimT')=p'{y), 4j{y,T')P, = p'{y), (3.13) 

where we set T' = (T — 1)/2 for the upper boundary hypercubes (recall that T has to be odd). 
The boundary four- component spinors p and p' are related to the boundary one- component 
spinors ( and (' as in eq. 



i i 

and analogously for p and p'. The action appropriate for the Schrodinger functional can now 
be written as ^ 

T'-l 

SsF = (2a)^E E Y.'^iy)i',T^My) + 4°^ + s'-P (s-is) 

y ?/4=i M 

with the boundary contributions 

3 

= (2a)^EE^(^'0)rfc^'^^(^'0) + (2«)'E^(^'0)^ir4Z^(^'l) 

y k=l y 

- (2a)3^^(y,0)f^^(y,0) (3.16) 
y 

and 

3 

S'^P = (2a)^^5^V^(y,r')f,I),^(y,T')-(2a)3^V^(y,r')Pof4^(y,T'-l) 

y fc=i y 

- {2afY.'^{y,T')Vli^{yX)- (3.17) 



y 



Using f| = Pif4Po — Pof 4P1 we see that the last term in both S'^'' and in involves a 
boundary field, and hence vanishes for homogeneous boundary conditions. 
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Miyazaki and Kikukawa already addressed the question as to whether there are addi- 
tional boundary counterterms contributing in the continuum limit. They would have to be 
operators of dimension three, i. e. of the form 

ASb = {2af J2 E mMy, 0) + cf ^m T')AMy, T')} . (3.18) 

y * 

Taking into account the discrete spatial rotational symmetry, parity — the projectors Po,i 
appearing in the boundary conditions (|3.13|) are invariant under parity — and the chiral 
U(l) symmetry of massless staggered fermions, Miyazaki and Kikukawa concluded that 

A, = f4, Ti (74®(747r) or (75®(757r), (3.19) 

where 7 = ^2^=113- Writing these possible boundary contributions in terms of the one- 
component fields X s-^id x they then argue that the last two terms contain derivatives and 
are of order 0{a) and therefore do not contribute in the continuum limit. Staggered fermions, 
however, have additional symmetries, overlooked in [§], namely shift invariance by one (fine) 
lattice spacing and a charge conjugation symmetry (see the appendix A for a summary of 
the symmetries of staggered fermions). Shift invariance (in spatial directions) excludes the 
last two possibilities in ( |3.19| ), making the more indirect argument in p] unnecessary. 



on the other hand, is excluded by the charge conjugation symmetry. Thus the only possible 
dimension three boundary counterterm is already present in the action. As mentioned above 
it involves a boundary field and it can therefore be absorbed into a renormalization of the 
boundary field, just as in the continuum (or Wilson fermion) Schrodinger functional 

3.1 Boundary counterterms at order a 



Though the bulk part of the action, eq. ( |3.15[ ), appears to have order 0{a) lattice effects 



— the aA^ part of — this is not so. The apparent 0{a) effect just stems from a "bad" 
choice in the construction of the four- component spinors in eq. (|3.4|) . It can be transformed 
away by using "improved" fields 

xliy) = xdy) - (^^^i.A^'^xdy)^ xliv) = xdv) ~ (^^^^.^xfiuiy)- (3.20) 

For our Dirac spinors this transformation becomes 

^\y) = Hy) - ■ ^-^(^)' ^'(^) = ^(^) - ■ ^i'^- (3.21) 

This choice of improved fields is not unique. We prefer a somewhat more symmetric form, 
which turns out to treat the fields near the boundaries in the Schrodinger functional more 
equally, 

V 

i,\y) = i;(y)-aJ2i^DM-l{Pt^-Pt^)- (3.22) 
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To apply the transformation, eq. ( |3.22| ), near the boundaries it is convenient to extend the 
fields %l) and to the region 7/4 < and 1/4 > T' by setting them to zero there. Inserting 
the improved fields into the action and doing a partial resummation one finds the additional 
term in the bulk [O 



-{2aYaY,^\y)DlVl,l^\y) 



which cancels the 0{a) term in eq. ( |3.15| ), up to higher orders in a. Luo [|T4| shows in addition 
that no dimension five operators (in the bulk) are allowed by the staggered symmetries and 
hence that no bulk 0{a) artefacts will appear at the quantum level. 

At the boundaries, however, 0{a) effects may occur. These come from dimension four 
operators at the boundary. Operators allowed by the staggered symmetries in the massless 
case under consideration are (see appendix A for a review of those symmetries) 



03 
C4 



fc=l 

3 



O2 = y^V^(7fc7475) ® {iblifT^ki^ 



k=l 

V^(0)f4- 

a 



:i)-^(o)]--[^(i)-^(o)] f4v^(o) 
:i)-^(o)] + ^[^(i)-^(o)] f^^(o) 



(3.23) 



(3.24) 



The term Oi and the combination O3 + already occur in the action, eq. ( p.l5[ ). O2 and 
the combination O3 — 04, both flavor symmetry breaking — they break the shift symmetry 
in the time direction which is broken already, of course, by the presence of the boundary in 
the Schrodinger functional — appear at the quantum level. 



3.2 Relation to the usual Schrodinger functional 

Starting from the transfer matrix for Wilson fermions, Sint arrived at fermionic boundary 
conditions that are different from eq. ( p.l3| ) For the four-flavor Dirac spinors considered 
here they would read 



P+V^(y,0) = p(y), ^(y,0)P_ = p(^), P_i;iyX)=p'{y), i^{yX)P+ = p' 



(3.25) 



with P4 



i ± r4 . Hence the Schrodinger functional with staggered fermions, deflned 
thus far, does not seem to agree with the "usual" deflnition of the Schrodinger functional 
in the presence of fermions. However, Sint noted that the staggered action and boundary 



^Yet other boundary conditions were considered by Symanzik |10|. However, those boundary conditions 
expHcitly break parity. 
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conditions in the classical continuum limit derived in [Q could be brought into the form he 
derived for Wilson fermions by a "chiral rotation" [|I3] 

i;' = R-ip and ^j' = ■ R (3.26) 

where ^ = ^4 ■ i?^ ■ ^4 and 

R = Rl{e^) = exp{t95i^l)} = cos^si - sin^gf^ = ^4(^5). (3.27) 



Defining 



one finds 



Po,im = R'M) ■ Po,i ■ {Rimr'. (3.28) 



^0,1(^^5) = I 



i±cos2^5(7475)® (7574)''Tsin2^5f4 . (3.29) 



Hence one can smoothly go from the "natural" Schrodinger boundary conditions for staggered 
fermions to the conventional ones for continuum (and Wilson) fermions, at least in the 
massless case under discussion. Since the rotation needed, eq. ( |3.27| ), is chiral this will not 
be true for massive fermions. We shall discuss massive fermions later on, in section ^. 

The "rotation" ( p.26| ) refiects an arbitrariness in the assignment of Dirac and fiavor 
indices in the construction of the four- component spinors in eq. (|3.4|). Such a fiavor-spinor 



rotation could be inserted directly into (|3.4| ). An allowed rotation should leave the kinetic 
term of the fermion action in the continuum limit unchanged. Thus we require that 

R'^ ■ ■ R-^ = V/i. (3.30) 

Therefore the generators of R need to be 

Ts(^Tp = [1, 75] ® [1, 7/., 75, i{l5l^,),'^{l^,lu)f■ (3.31) 

Note that all these generators are hermitian. In particular the generator iT^, and hence the 
rotation Rl, ( p. 27] ), is allowed. 



4 The Schrodinger couphng in the massless case 

The first real application of the Schrodinger functional formalism in lattice gauge theory 
was the computation of a well-defined renormahzed coupling The special boundary 

conditions allow introduction of an external field to probe the system, and the finite size of 
the system gives a definite scale. In this section we will compute the one-loop contribution 
from the staggered fermions to this coupling and show that it agrees with the result using 
(improved) Wilson fermions or a continuum regularization. This serves as a test for the 
correctness of the Schrodinger functional setup with staggered fermions. As a by-product 
we will also obtain the contribution from the staggered fermions to the pure gauge 0{a) 
boundary counterterm with coefficient q in the notation of |]T],§|,|§|. 
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The external field is introduced via the boundary gauge fields taken as abelian with 



Wk{x) = diag (e 



i<t>i/L „i4>2/L i4>-j./L 



Wl{x) = diag(^e''^'i/^,e*^2A,e''^3A 
They lead to the classical gauge fields 



0, 



(4.1) 



(4.2) 



As [0,0 we choose the boundary fields to depend on a parameter, which we denote by u, 
through 



TT 



3 

-TT — UJ, 



1 vr 1 

2^' '^^=3-2"' 



TT 1 
1 CO, 

3 2' 



27r 1 

h -UJ. 

3 2 



(4.3) 



The "Schrodinger functional coupling constant" is then defined as 

2 



k=l2\- 



w=0 



sm + sm —-= 

\3LT J \^LT ) 



(4.4) 



with T = L such that (7^ depends only on one scale, = g'^{L). The normalization k has 
been chosen such that g equals the bare coupling at tree-level without any cutoff effects. 

The one-loop contribution from the staggered fermions to the coupling constant eq. (O) 
comes from the derivative of the fermion fiuctuation determinant. The fermion boundary 
fields here are set to zero. Then the fermion action can be written schematically as 



T-l 



(4.5) 



X X4 = l z 

As usual with staggered fermions, one easily sees that 

(-l)l-lM.,,(-l)N = Mt . 



(4.6) 



Thus (— 1)'^'Mx.,2 = Aix,z is hermitian, and has the same determinant as M^^z- The one-loop 
contribution to the running coupling is given by [S] 



9^ = 9o+ Pi9o + C'(^o)' Pi = Pho + nfPi,i, 
with the fermionic contribution 

1 d 



Phi 



log det M 



w=0 



(4.7) 



(4.8) 



Here rij = 4 for one flavor of staggered fermions, since they correspond to four flavors of 
continuum fermions. 
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The eigenfunctions of A4 are of the form 



X{x) = e^(^^+"")Va(x4), Ux^ = 0) = /5(X4 = T) = 0, (4.9) 
where a is a 3- dimensional vector with = 0, 1 and 

Pfc = -^ + -^, n, = 0,... ,-Lfc-l. (4.10) 

As usual for staggered fermions the momentum components go only over half the Brillouin 
zone interval. The remainder, afcvr, becomes an "internal" index, loosely corresponding to 
spin/flavor. 

Introducing r]^^^ as 

/ N I 1 for z/ < u , , 

= n f > (^-^^^ 
10 tor z/ > /i 

such that the staggered phase factors become r]^(x) = (— l)'?'^'''^ and realizing that (— l)'^' = 
(^_l-jv''*^x+x4 gj^^ fdi^i) satisfies 

(-l)^^-^«,/f/|^(^4) = lfl\x, + 1) - X^^sin(r« + ^^^^x,)i-ir%^^^,,^^^^,jf^^\x,) 

k=l 

- ^/|Ha^4-l). (4.12) 

Here 6 means 6 modulo 2, r^*'' = Pfc + and = (0^ — (j)i)/{LT). 

Essentially, due to the species doubling of staggered fermions in the time direction, even 
though we are dealing with fermions, eq. ( |4.12| ) is a (hermitian) second difference equation. 



We can therefore directly apply the recursive method of appendix B of to compute the 
fluctuation determinant for fixed p. 

To be precise, the coupling is defined in as eq. ([4.4|) with T = L so that it depends 



on a single scale, L. Unfortunately this is not possible for staggered fermions, since, as we 
have seen, L/a must be even but T/a must be odd (we have restored here the dimensions 
of L and T).^ The closest we can get to T = L is therefore T = L ± a. They all coincide 
in the continuum limit, but at finite a taking T = L ± a introduces additional 0{a) effects. 
Averaging the couplings obtained with T = L + a and T = L — a cancels this additional 
0{a) effect. 

We have evaluated p^{^i of eq. ( [4.8|) for L/a ranging from 4 to 64, in steps of 2. Here 
the superscript ± stands for the choices T = L ± a. One expects pi^i{L/a) (we omit the 
superscript for the generic case) to be given by an asymptotic series of the form 

Pi,i(L/a) = ro + So log(L/a) + (n + si \og{L/ a)){a/ L) + 0{\og{L/ a){a/ Lf). (4.13) 



•^The fermion fluctuation determinant would be zero for T/a even, yet another confirmation that T/a 
must be odd. 
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The first few coefficients in eq. ( [4.1^ ) can be extracted by first cancelling higher order terms 
in a/L through numerical differentiation and then checking for stability as L/a is increased 
0. So, the coefficient of the logarithmically divergent term in the continuum limit, should 
just be 26o,i = — l/(127r^), the fermionic contribution, per flavor, to the /9-function 0, and 
thus absorbed by renormalization. We indeed found this result to about 6 digits accuracy 
for all cases considered. 

si was found to be compatible with zero, albeit only with an accuracy of about 4 digits. 
To extract tq we assume the exact value for sq; to obtain ri we assume in addition that 
Si = 0. Then we found the values listed in Table ^ Tq obtained from the choices T = L ±a 
always agreed to the accuracy given. 








J+) 
' 1 


' 1 


' 1 





-0.004416(1) 


-0.02003(1) 


0.03897(1) 


0.00947(1) 


7r/5 


-0.00579695(2) 


-0.023113(5) 


0.042061(5) 


0.009477(5) 


1.0 


-0.0068642(1) 


-0.02552(1) 


0.04447(1) 


0.00947(1) 



Table 1: The first two 'non-log' terms in the expansion eq. ( [4. 13|) of pi^i. comes from 
the choices T = L ± a and rf"^^ is the average of the two. 

To of eq. (|4.13|) is the finite part of the fermionic contribution to the one-loop relation 
between bare lattice and running Schrodinger functional coupling. Using the known one-loop 
relation between lattice and MS coupling we can obtain a relation between the Schrodinger 
functional coupling a and a^^g 

c^mil) = + ci«'(g) + 0(«'), a{q = 1/L) = f{L)/{An). (4.14) 

This relation between two renormalized continuum couplings has to be independent of 
the regularization used to obtain it. Thus, we should find the same result for ci^i in 
ci = ci,o + n/ci,i as [g. ci,i is given by 

Ci,i = -47r[P4 + ro] (4.15) 

where P4 = 0.0026247371 is the finite contribution from staggered fermions to the one-loop 
relation between lattice and MS coupling, computed in [|ll|.Q Using tq from Table |I] we 
obtain 

fo.02251(l) for 9 = , , 

cii = < ^ ' 4.16 

[0.0398632(2) for = 7r/5 

in good agreement with the results of Sint and Sommer for Wilson fermions |^ . This confirms 
that our Schrodinger functional for massless staggered fermions is a correct regularization of 
the continuum Schrodinger functional with massless fermions. 

''The better accuracy quoted here for P4 is obtamcd from eq. (6.12) of Q using the more accurate values 
for Pi and P2 from . 
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4.1 Lattice artefacts 



The vanishing of si indicates the absence of bulk 0{a) artefacts 0, as expected for staggered 



fermions [Q. The non- vanishing of ri reflects the presence of boundary 0{a) effects. As 
mentioned before, the fact that we are not allowed to take T = L for the computation of the 
running coupling introduces additional 0{a) effects that can not be cancelled by boundary 
counterterms. They can, however, be cancelled by averaging over the choices T = L ± a, 
as can be seen from Table |l]:0 r(""^ is independent of the spatial boundary conditions 



parameterized by 6. This remaining 0{a) effect can be absorbed into the pure gauge 
boundary counterterm (q — 1) X]fc=i ^^(-^fc4-^fe4) [§] by choosing q = 1 + (cj^''^'* + nfcf''^'')g'^ 
with the pure gauge part c[^'°'* = —0.08900(5) and the fermionic contribution 



.(1-1) _ 



0.00474(1). (4.17) 



After cancelling the 0{a) boundary lattice artefact by the boundary counterterm, higher 
order lattice artefacts, both from the bulk and the boundary, remain. With our results we 
can study them for the fermionic contribution to the step scaling function [^,3, S(s, u,a/L) 
which is the coupling g'^{sL) at scale sL when keeping the coupling at scale L fixed at 
g'^{L) = u. In the continuum, the one-loop fermionic contribution to a{s,u) = S(s,m,0) is 
26o,i^/ log(s) u^, while on the lattice it is given by nf[pi^i{2L/a) —pi^i{L/a)] u^. We take 
s = 2 and compare the lattice result (per continuum flavor) — we consider here the average 
between the choices T = L ± a — with its continuum limit: 

p[-\2L/a)-p[fiL/a) 

= — 2k;:j^2 — • ('-^'^ 

Deviation of 6i^i{a/L) from 1 at flnite a/L is a lattice artefact. Note that we deflne the lattice 
artefact from the fermions, 6i^i{a/L) — 1, with respect only to the fermionic contribution in 
the continuum limit, in contrast to 6i^i{a/L) is shown in Figs. |^ and for 6 = and 
vr/S, respectively. Shown in both cases is the result with and without the cancellation of 
the 0{a) part of the lattice artefact by the pure gauge boundary counterterm. As can be 
seen from the two cases, the higher order lattice artefacts can depend very sensitively on the 
observable considered, here the step /3-function for different values of the spatial boundary 
conditions. After cancellation of the 0{a) part by the boundary counterterm, one expect the 
remainder to go asymptotically (up to logs) like (a/L)^. Rough estimates of the coefficient 
of the {a/LY contribution from the flgures are -7.5 and -0.5 respectively. For 9 = 1 the 
coeflicient becomes about 4.5. These rather large variations should serve as a caution to 
drawing conclusions about the order of magnitude of cut-off effects from studying just one 
observable. 



^We expect this to be true also for the pure gauge part that, in a simulation with dynamical staggered 
fermions, would obviously have to come from an average of systems with T = L ± a. We have not checked 
this explicitly. But S. Sint has checked that the statement holds for the one- loop contribution from Wilson 
fermions. In addition we have performed a pure gauge MC simulation for /? = 5.9044 with L/a = A and 
T = LzLa and verified that the average of the inverse of the non-perturbatively computed coupling constants 
agrees, within errors, with the result for T/a ~ L/a = 4 in ||]. 
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B = 0.0 




0.00 0.01 0.02 0.03 

(a/L)^ 



Figure 1: The ratio of the fermionic contribution to the step scahng function at one-loop 
from the lattice to its continuum limit, ( |4.18| ), for 6* = 0. Crosses show the result obtained 
from the fermion fluctuation determinant without cancellation of the 0{a) part by the pure 
gauge boundary counterterm c^t'^\ and octagons the result after the cancellation. The line 
shows the 0{a) part of the lattice artefact that is cancelled by the pure gauge boundary 
counterterm. 



5 The Schrodinger functional for massive staggered fermions 

We have seen that with a chiral rotation, eq. ( p.27[ ) with 9^ = 7r/4, we can turn the "natural" 
staggered boundary conditions, in terms of four-component Dirac spinors, into the usual 
boundary conditions obtained when starting with Wilson fermions. Since this is a chiral 
rotation, the mass term, omitted so far, will not remain invariant. Thus, to obtain the 
usual Schrodinger functional boundary conditions with the usual mass term we have two 
choices. Either we change the boundary conditions for the staggered fermions, or we use an 
unconventional mass term. 

Here we consider the second alternative. Thus, instead of the usual mass term mipip in 
terms of the four-component spinors before the chiral rotation we consider rm/jT^ip, which 
becomes the usual mass term after the chiral rotation. Using eq. (p.4|) this mass term can 
be written in terms of the one-component spinors as: 

Sm, =mY,V4{xmx) 1^ [1 - (-1)^1 Ulix - l)x{x - 4) - 1 [1 + (-1)^*] U,{x)x{x + 4)| , 

(5.1) 
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1.2 



= tt/5 




0.9 I — L 
0.00 



0.01 0.02 0.03 

(a/L)^ 



Figure 2: Same as Fig. |1|, but for 6' = 7r/5. 



where we included the gauge fields to make this non-local (on the fine lattice) mass term 
gauge invariant. 

However, this new mass term, eq. ( |5.1| ), introduces an 0{a) term in the bulk. This can 
be seen by writing it in terms of the improved four-component spinors, eq. (|3.22| ), 



To the same order, we can replace D4 in the second term by V4. Subtracting this term, we 
arrive at the 0{a) improved mass term 



Yl ^(-l)"'^4(x)x(x) (U,{x)x{x + 4) + Ulix - l)x{x - 4)) . (5.3) 



A few comments are in order here: the mass term ( |5.3|) obviously singles out the time 
direction. But the time direction is singled out in the Schrodinger functional formalism 
anyway. This mass term is invariant under the "chiral" U{l)e symmetry, but it breaks the 
shift symmetry in the time direction. Indeed, a shift by one lattice unit in the time direction 
reverses the sign of the mass. Hence, for the usual setting with (anti-) periodic boundary 
conditions, the model is invariant under a change of sign of the mass. The free spectrum, for 
example, is a function of m^. However, such a shift in the time direction is not possible at 
the boundary. The Schrodinger functional is thus sensitive to the sign of the mass, and we 
therefore expect the Schrodinger coupling to depend linearly on the mass for small masses. In 
the continuum limit, the shift symmetries become part of the spin and vector and axial fiavor 
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symmetries. Thus breaking of the shift symmetry imphes breaking of the chiral symmetry 
in the continuum, just as is expected of a mass term. 

Since a mass term can break a symmetry, counterterms have to be re-examined in its 
presence. The mass term ( |5.3| ) only breaks the shift symmetry in the time direction, but 
leaves the "chiral" symmetry unbroken. In the bulk, the only allowed counterterm 

multiplicatively renormalizes the mass. On the boundary, the only 0{a) term allowed in 
addition to those in eq. (|3.24|) is 



(5.4) 



i.e., the new term now already appearing in the action. 

Having found a mass term that is designed to correspond to the usual mass term in 
the Schrodinger functional formalism, we can check this at one-loop by computing the 
contribution to the running coupling. The computation goes as in the previous section. 
The eigenfunctions now satisfy 



3 

- 5^^sin(r« +#X4)(-l)"^5„.^^^4)+,^.),^^/f (X4) (5.5) 



fc=i 
1 



- ^[^-i-irrn]fjl\x,-l). 

Since the scale for the running coupling is given by the spatial system size L we want 
to keep z = mL fixed when varying L/a. Again, we compute the contribution pi i{z, L/a) 
for L/a = 4 up to 64 and extract the leading terms in the asymptotic expansion eq. ( [4.13| ). 
A few examples are listed in Table ^j. Again we find sq = — l/(127r^) to our accuracy, and si 
compatible with zero.0 



e 


z 


^0 


J+) 
' 1 


J~) 
' 1 


' 1 





1.0 


0.0024408(1) 


-0.00426(1) 


0.02319(1) 


0.00947(1) 


7r/5 


1.0 


0.00179605(3) 


-0.00617(1) 


0.02511(1) 


0.00947(1) 





2.0 


0.0064310(1) 


0.00223(1) 


0.01671(1) 


0.00947(1) 


7r/5 


2.0 


0.00632882(5) 


0.00203(1) 


0.01691(1) 


0.00947(1) 



Table 2: tq and ri's as in Table |T] but in the massive case with z = mL held fixed. 



The difference 01,1(2;) — ci,i(0) = — 47r[ro(-2)— ro(0)] should be regularization independent. 
Comparing with the Wilson fermion results of we indeed find agreement within the 
accuracy given in Tables |l] and |^. 

^We actually discovered the presence of a bulk 0{a) contribution for the unimproved mass term, eq. (^.l|), 
by finding a non- vanishing si. 
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The value of r} is independent of z and hence cancelled by the pure gauge boundary 
counterterm with coefficient ( ^4.17|) . No mass dependent 0[a) bulk term of the form 
m\x{F^vF^y) is needed, in contrast to the case of Wilson fermions 0. 

If we do not insist on reproducing the conventional [i.e. Wilson fermion inspired) 
massive Schrodinger functional in the continuum limit, we can use the "natural" staggered 
boundary conditions and the usual (degenerate) staggered mass termQ without chiral rotation 

S^ = rnY^ X{^)x{'^) = (2a)'^m ^'ip{y)ip{y). (5.6) 

X y 

This leads to a well-defined Schrodinger functional in the continuum, albeit with boundary 
conditions that mix the four continuum fiavors. In contrast to the "conventional" case now 
both action and boundary conditions are invariant under m — > —m and thus the Schrodinger 
functional running coupling defined by these conventions is an even function oi z = mL. 
Inclusion of the mass term ( |5.6|) into eq. ([4.12|) is straightforward, giving 

i-^rM^^^,ff{x,) = \ffix, + 1) - X^.sin(r« + V^«a:4)(-l)"'=5,^^.,^4)+,^(.,,,^/f (X4) 

fc=i 

+ "^W),/?/|^(^4)-^/l^(^4-l). (5.7) 

Proceeding as before, we now obtain the sample results listed in Table |^. We note in 
particular that r^"^"* is unchanged from before and hence cancelled by the pure gauge 
boundary counterterm with coefficient ( [4.17|) . 
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z^ 


^0 


.(+) 
' 1 


' 1 


' 1 





1.0 


-0.0004189(1) 


-0.00853(1) 


0.02747(1) 


0.00947(1) 


7r/5 


1.0 


-0.00195976(3) 


-0.013405(3) 


0.032350(3) 


0.009473(3) 



Table 3: tq and ri's as in Table |I]but in the massive case with the standard staggered mass 
term, (|5^). 

The mass term (|5.6| ) breaks the "chiral" U{l)e symmetry. As a consequence the 0{a) 
boundary counterterms 

Oy = m^/^(7475) ® (7574)^^^- (5.8) 

are allowed in addition to those in eq. ( p.24| ) and the term in ( |5.4| ). 

^We do not consider mass terms that break some of the staggered symmetries (other than the chiral 
symmetry) which might be used to hft the degeneracy among the four continuum flavors. 
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6 Conclusions 



We have discussed the Schrodinger functional with staggered fermions, extending and 
completing previous work by Miyazaki and Kikukawa. In particular, we have shown that 
for massless fermions the Schrodinger functional constructed agrees, in the continuum limit, 
with the conventional Schrodinger functional for fermions as described by Sint, by examining 
the Schrodinger functional for staggered fermions in terms of four flavors of Dirac spinors. 
We found one boundary counterterm that contributes in the continuum limit, just as for 
Wilson fermions and for continuum fermions with dimensional regularization. It can be 
absorbed by a renormalization of the boundary fields. 

We have computed the one-loop contribution from massless staggered fermions to the 
running Schrodinger functional coupling and found agreement, in the continuum limit, with 
the results of Sint and Sommer for Wilson fermions. We computed the one-loop contribution 
from staggered fermions to the pure gauge 0{a) boundary counterterm, which cancels the 
0{a) lattice artefact in the running coupling. We discussed the higher order lattice artefacts 
in the step scaling function at one loop. They can be as large as 20% on an L/a = 6 lattice, 
and they can depend rather sensitively on the observable considered, in our case the coupling 
defined with different fermionic boundary conditions. 

A chiral, flavor changing rotation was needed to bring the boundary condition of 
staggered fermions, expressed in terms of four flavors of Dirac fermions, to the conventional 
form. The mass term is not invariant under such a rotation. We therefore considered an 
unconventional mass term for the staggered fermions, constructed to give the conventional 
Schrodinger functional for massive fermions in the continuum limit. The construction was 
verified by computing the one-loop contribution to the Schrodinger coupling constant for 
fixed z = mL and reproducing results by Sint and Sommer. 

One unpleasantness in defining the Schrodinger functional running coupling constant 
with dynamical staggered fermions is the fact that we cannot take T = L, since L/a has 
to be even while T/a has to be odd. Therefore the coupling does not strictly depend on 
a single scale. Taking T = L ± a gives a single scale in the continuum limit (a — ^ 0), but 
introduces 0{a) effects at finite lattice spacing. As discussed, these 0{a) effects, which can 
not be cancelled by boundary counterterms, can be cancelled by averaging over the coupling 
obtained with T = L ± a. However, since simulations with staggered fermions tend to be 
less costly than simulations with Wilson-type fermions — staggered fermions have four times 
fewer degrees of freedom, and no fine-tuning is required to make them massless — having 
to do two simulations for each L/a might not be such a big price to pay. 

Having established the Schrodinger functional with dynamical staggered quarks opens 
the possibility to use it for the non-perturbative computation of improvement coefficients 
and current renormalization constants for gauge field ensembles with dynamical staggered 
fermions, just as it is done for quenched simulations by the ALPHA collaboration We 
envision here the computation of the clover coefficient and current renormalization constants 
for (improved) Wilson valence quarks. These can then be used to compute phenomenolog- 
ically interesting quantities, such as e.g fs and /d, with Wilson valence fermions on gauge 
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configurations generated with dynamical staggered fermions |]TB[ . Comparison with quenched 
results then allows an estimation of quenching errors in those quenched computations. 
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Appendix A 



In this appendix we briefly review the symmetries of staggered fermions |2y|T^. Unless 
otherwise noted, the one-component spinors transform under a symmetry transformation T 
as 



X{x) x{.x) = r]TiT'^x)xiT~^x). 



The gauge fields transform as 



U{x,z) U{T-'x,T-'z) 



(A.l) 



(A.2) 



with 



U (x, z) 



U^{x) for 2; = X + /i 
f/^(2;) for X = 2; + /i. 



(A.3) 



(i) Rotations by 7r/2: 

Here we consider rotations around the center of a hypercube (p < a): 

Xa -Xp + l (A.4) 

Xf, Xf, for fij^p,a. 

Then we have 

In terms of the hypercube block coordinates: 

Vp VcT, Va -> -Vp, Vp -> Vp for fi^ p,a, (A.6) 
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and the fields ip and ip transform as 

^{y) ^ \{'^ + ipi'r)®{il-i'^p)-^{R-Hy) 

i,{y) ^(i?^i^).i(l-^^^,)®(^J-^J). (A.7) 

(ii) Reflections on hypcrplancs: 

Here we consider reflections on hyperplanes through the center of a hypercube: 

Iff • Xp > Xp -\- 1 



X 



Xp, for // 7^ p. (A.8) 



Then we have 

rji,{x){-ir^+-+-\ (A.9) 
In terms of the hypercube block coordinates: 

Vp^-Vp, Vi^^Vn ioridy^p, (A.IO) 
and the fields ip and ip transform as 

^(y)^{ipi5)®il -^iH^y), ^(y)^^(iH^y)-(i5ip)®iI- (A.ii) 

Combining 3 reflections, one orthogonal to each spatial direction gives the parity 
transformation 

P : y^-y, y^ ^ y^, (A. 12) 

with 

i^{y) ^ 74^75 -i^iPy), i^{y)^i^{Py)-i4®iI- (a.13) 

(iii) Shift invariance: 



Xp Xp for A* 7^ p. (A. 14) 



Then we have 



ry^(a;)(-l)^''+^+■••+^^ (A.15) 

and the fields ip and ip transform as 

1 1 

V'(y) ^ ^Sy,y'{i ® tJ - 7p75 ® tD • '^{y') + 2^y+P,y'('^ ® tJ + 7p75 ® ll) ■ ^{y') 

i^iy) ^ lSy,y4iy') • (1 «) 7j + 7p75 ® ll) + \5y,y'~pW) " (1 7J - 7p75 ® lI)- (A.16) 
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Symmetries (i) through (iii) are the space-time symmetries of staggered fermions. In 
addition we have the invariances: 

(iv) U(l)-invariance: 

Xix) ^ e-x(^), X{x) ^ e-'"x{x), (A.17) 

which just becomes 

^(^)^e^"^(y), v5(y)^e-^"V5(y)- (A.18) 

and 

(v) U(l)e-invariance: 

^ e'^<'\{x), x{x) - e'^<^^x{x), (A.19) 
where e{x) — (—1)1^1 This becomes the "chiral" transformation 

^{y) e^P'^^®^^ . ip{y), ^{y) ^ tp{y) ■ e^^Ts^Ts"^. (A.20) 

This chiral symmetry protects the zero-mass hmit for staggered fermions, since the usual 
mass term mJ2xX{x)x{x) = m2'^J2yi'{y)'^{y) is not invariant. 
Finally, staggered fermions have the discrete 

(vi) Interchange symmetry: 

X{x) e{x)x^{x), x{x) -x^{x)e{x), (A.21) 

where T stands for transpose (as a color vector). The gauge fields transform as 

U,{x) ^ U;(x). (A.22) 

In terms of the four- component spinors this becomes charge conjugation symmetry 

^ (A.23) 

Here C = C ® (C^^)^ where C is the usual Euclidean charge conjugation symmetry matrix 
satisfying 

Cl,C~' = -7j 

Ci^l^C-^ = (757^)^ (A.24) 

-C = C'^ = C-^^Cl (A.25) 
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Appendix B 

In this appendix we sketch the construction of a "two-time-shce" transfer matrix for 
staggered fermions, following quite closely the construction in |TT| for the "reduced" 
staggered fermions. We consider staggered fermions with the standard (flavor degenerate) 
mass term. The partition function is 

Z = J [DU][dxdx] exp{-SG - Sp}. (B.l) 

with Sg the usual Wilson gauge action and 

Sf = ^Y1 [U^,ix)x{x + Ul{x - fi)x{x - f^)] mx{x)x{x). (B.2) 



We now change integration variables for each pair of time slices X4 = 2r and 2r + 1 

al{x) = Pe{x)^r]4{x)x{x,2T) + Po{x)x{x,2t) 

ar{x) = Pe{x)x{x,2r+l) + Po{x)^vdx)x{x,2T + l) 

(3l{x) = P,{x)x{x,2T) + Po{x)^r],{x)x{x,2T) (B.3) 

Prix) = Pe{x)^V4ix)x{x,2T + l)+Poix)xix,2T + l) 

where Pe,o{x) = | (l ± (— l)'"^') are the projectors onto even and odd spatial sites x. We 
have used here that staggered phase factors r/^ are independent of the time coordinate; in 
particular, ?74(x, 0:4) = ?74(x) = (— l)'^'. 
Introducing 

fjkix) = r]i{x)r]k{x), fjk{x + k) = -fjk{x), 
we can write the staggered action in temporal gauge, Ui{x) = 1, as 

T X 

+ P,{x)J2[4i^)US)Uk{x,2T)aUx + k)+Pl{x + k^^^^ 

k 

+ar{x + k)fjk{x)Ul{x, 2t + l)ar{x) + pr{x)fik{x)Uk{x, 2t + l)pr{x + k) 
+ Po{x) J2 [pKx)f]k{x)Uk{x, 2t)PI{x + k) + alix + k)fik{x)Ul{x, 2r)a\{x) 

k 

+(3r{x + k)fjk{x)Ul{x, 2t + l)Prix) + ar{x)rikix)Ukix, 2t + l)ar{x + k) 
+ 2m[a\{x)(3l{x) + (3r{x)ar{x)\] . (B.4) 
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Next we introduce operators a{x), a^{x), P{x) and with the anticommutation 

relations 



(B.5) 



and all other anticommutators vanishing. On the Fock space, spanned by these operators, 
we consider the coherent states 



ttr, /3r^ = exp < |^d^(x)Q;^(x) + (3^{x)Pr{x) > |0) 



(B.6) 



and 



al, Pi = (0| exp i J2 \(yl{x)a{x) + Pl{x)(3{x) 



(B.7) 



where |0) is the Fock vacuum. It can be shown that these coherent states satisfy the 
completeness relation 

l = j [rf^^a,] [rf/?^d/3,]exp|-^ [4(^)ct,(^)+^t(^)^^(^)]| ar,Pr)(^al,Pl , (B.8) 

and that for a normal ordered operator A =: A{d'^ , a, /3) :, where normal ordering means 
that all daggered operators appear to the left of all non-daggered operators. 



X) 



(B.9) 

We can see that the terms in the exponentials in ( p.8|) and ( |B.9| ) reproduce the terms 
from the first line of the action 5*^, eq. (p.4| ). The other terms can be reproduced from the 
transfer matrix (in temporal gauge) 



:b.io) 



with Tg the pure gauge transfer matrix (in temporal gauge) ||21|| — we need a total of 
since we move two time slices forward. For Ti;' we find 



exp 



E(^. 



[X 



a{x + k)r]k{x)Ul{x)a{x) + P{x)rik{x)Uk{x)/3{x + k) 

+Po{x) \f3{x + k)fik{x)Ul{x)p{x) + a{x)fjk{x)Ukix)a{x + k) 
+m(3{x)a{x 



[BID 



With some algebra one can show that TtXi;' is equal to the square of the "one-time-slice" 



transfer matrix of |TT[ when gauge fields and the mass term are included in the latter. 
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The "two-time-slice" transfer matrix T, ( [B.1U| ), is self-adjoint and positive. Therefore, 
restoring the lattice spacing in the time direction, a^, we can define a Hamiltonian by 



T = e" 



-2atH 



(B.12) 



H is a complicated function of the fermion and gauge field operators. It simplifies in the 
naive time continuum limit, at 0, where we can neglect contributions of 0{af) and higher. 
The fermionic part becomes in this limit 

Hp = '^l^Pe{x) a\x)fjk{x)Uk{x)d\x + k) + a{x + k)fik{x)U^^ 

X 

+I3\x + k)fik{x)Ul{x)f3\x) + (3{x)fikix)Ukix)(3{x + k) 

+ ^Po{x) [p^{x)US)Uk{x)P^{x + k)+${x + k)U^^^^^ (B.13) 

+a\x + k)'fik{x)JJl{x)a\x) + a{x)'fik{x)Uk{x)a{x + k) 
a\x)j3\x) + j3{x)a{x) 



+ 



We can bring this Hamiltonian to a more conventional form with the canonical transforma- 
tion 

= -^Pe(f) - a(x)) + -^Po(x) - 



(B.14) 



The operators x ^i-nd ip and their adjoint's and ip'^ satisfy anticommutation relations like 



eq. (|B.5|) . Expressed in terms of these operators the Hamiltonian (|B.14|) reads 
Hp = ^ ^jj'^ {x)fjk{x)Uk{x)ij^ {x + k)+ ^{x + k)fjk{x)Ul{x)ij{x) 

x,k 

+x{x)fik{x)Uk{x)x{x + k) + x\x + k)fik{x)Ul{x)x\x) (B.15) 
+ m^(-l)'^' ■ip\x)-ii){x) - x\x)x{x) 



We recognize ( |B.16| ) as the Susskind Hamiltonian |^2[ for two independent staggered fermion 
fields ip and x- 
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